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AND SHOCK CONDITIONS FOR GASES* 
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This paper is concerned primarily with the mathematical structure 
of the derivation of the fundamental differential equations of fluid 






dynamics and the conditions for shock waves in gases. No attempt is made 
to treat the existence theoretic aspects of the problem or to discuss 






underlying physical ideas. 







VISCOUS AND COMPRESSIBLE FLUIDS 










Fundamental dynamical assumptions. A fluid (liquid or gas) wil] 


he considered from our present standpoint as forming a continuous 






medium. Let the motion of the fluid be represented by equations of the 





a 


; 1 a 7 s 
form x = f (x,t) where the x and x are coordinates of the same 





| rectangular system and t denotes the time. A particle, initially at the 






point x, will be located at the point x at time t in accordance with 





a 1 . . a 
=» these relations. It will he assumed that the equations x = f° (x, t) have 


4 unique inverse at any time t and it suffices for our purpose to 





a - . . 
suppose that the functions f (x,t) are continuous and have continuous 





partial derivatives through the third order with respect to the initial 
a 5 ° a 
coordinates x and the time t. The velocity components u (x,t) are then 






ce: » : . a . Ase) 
lefined in the usual manner by the equations u If" /9t and these 





2 = . . a “ 
quantities will be continuous in x and t and will have continuous 






® derivatives through the second order with respect to these variables. 
Now consider a finite volume V of the fluid and denote by S the 





boundary or surface of V. We suppose that V is a continuous (1,1) may 







e of a cube and that S is a closed regular surface ; these conditions wil] 
» be retained under motions of the fluid particles by which V and S are 
» composed in consequence of the above assumptions on the functions f 

© All volumes and their boundaries in the following discussion are assumed 
= implicitly to be of this type. We make the following dynamical assump- 





= tions. 
(a) The rate of change of the momentum of V parallel to any fixed 
direction is equal to the component parallel to this direction of the 






total external forces acting on V. 
(6) The rate of change of the angular momentum of V about any fixed 
line is equal to the moment of the external forces acting on V about 
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this line. be 
We assume also the law of the conservation of mass which states that de 
the mass of the moving volume V remains constant. This condition can be 
expressed by writing (4 
(1) 5 {edv = 0, ™ 
' ou 
where ( is the density or mass per unit volume at points of the fluid. 2 
In the case of compressible fluids for which ( is not constant the 
above assumptions must be supplemented by an equation of state; this 
consists usually of a relation between the density and pressure (defined ti 
in § 6). 
The assumptions (a) and (/) are fully satisfied if the conditions 
involved are expressed for each of three fixed and mutually perpen- 
dicular directions in space. We take these perpendicular directions @: 
to be the axes of the rectangular coordinnte system x to which the ” 
position of the fluid particles is referred. Relative to these coor- 
dinates let X°* be the components of an external or applied force per (5 
unit mass acting in V. Also let T be the components of a surface force 
per unit area over the surface S. The momentum per unit volume in the Eqn 
fluid will be given by ~u*, Finally we need an expression for the moment th 
of a force F about the coordinate axes. But the moment of a force F 
about the x° axis is given by M, = Cony x5F” where F” are the components fF eg 
of F and the quantities e,,, are defined by the following two require- 
ments, first e,2, = 1 and second the e,,, are skew-symmetric. Using 
these designations, conditions (a) and OD have their formal expression 
in the following two sets of equations 
(2) 4 [ pu*dv = [ Tus + { oxy, 
dt -y s v 
(3) 4 J apy 2u%dV = [ egg, xT'dS + [ egg, xX 70d. os 
V S V 
Th 
In the following it will be supposed that X = 0, i.e., no applied 
forces act on the fluid. Also it is assumed, preceeding $18, that all (6 
functions of the variables x* and t which enter into our discussion are 
continuous and have continuous partial derivatives for a!! values of Th 
these variables under consideration. Finally, we remark that covariant by 
and contravariant indices have been interchanged at will as is permiss- 
able since we deal exclusively with rectangular coordinate systems. 
Repeated indices in any term are to be summed in accordance with the th 
usual convention. un 
2. Differentiation of integrals with variable boundaries. Let V(t) (7 
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be a volume in the x space which may vary with the time t and let S 
denote the surface of V. It can be shown that 







Sf Peay [ faves GdS 
(4) hey I ~ Sy at — 






where G is the component of the velocity of the surface S along the 
outward normal to S. In the particular case that the variation in V is 
determined by the motion of the fluid particles we have G = ug’* where 
y* are the components of the unit normal vector directed outward from S. 







3. Equation of continuity. Taking f = 9 and G= u,l™ the applica- 
tion of (4) to the condition (1) yields, 





0p 
9? + Puav*dS = [ ay + (opu.) ,dV= 0. 
|= J, i vot y oe 





Since V is arbitrary it follows that, 








0p 7 dp 
(5) +t + (Pua) a - rn + Puasa = 0. 







Fquation (5) is called the equation of continuity; it is equivalent to 
the condition (1) for arbitrary volumes V. 

A useful result can be derived from (4) when G = u,v® and the 
equation of continuity is satisfied. We have 











J fugPdS = Si (fuga = | fauad¥ + J fuaadv = 






f do 
fhanav-[F% g dV. 
df 


d 
Hence, — J fdV = V - [f = dV, 
at? Vdt Vodt 






Then, replacing f by Pf, we obtain, 


d df 
(6) — dV ef — dV. 
dt J, Pf vy” at 





This equation is satisfied in a fluid when the variation in V is produced 
by the motion of the fluid particles. 












4. The stress tensor. It can readily be shown’ in consequence of 
the condition (2) that the above vector T representing the force per 


unit area or stress over the surface S has components of the form 







(7) Ty = ygvF, 
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where the 7,g are functions of position and time. From the vector 
character of T and v and the relation (7) we infer that the quantities 
tg are the components of a tensor under coordinate transformations. It 
will now follow from (2) and (3) that the tensor 7 is symmetric. 

Using (6) equations (2) and (3) become 


p — dV = [ tevfdS = [ 7 av, 
Jy dt ‘s °F per 


ee | 
Li qe eal > { egpyxht,s v?dS = [, (eapyx ys ),s dV, 


 . 


When we carry out the diflerentiations indicated in the second of these 
equations the two sets of equations can be written 


(8) | cae 
lOoO— -f dv = 0, 
Ny dt aB,B 
6) [ B du* ; . 
(9) ty CaBy* (0 de _ Ty8 5) dV + y Ca By By dl = Q. 


Then, since V is arbitrary, we have 


; du* 
10) PS * eB S 
from (8). But from (9) and (10) it now follows that CanyT py = 0 and 


these conditions imply the symmetry of the quantities 7,2. Hence we see 
that (2) and (3) can be replaced by the equations (10) with the condi- 
tion that the tensor T be symmetric. 


5. The distortion tensor. Suppose that in a smal?! interval of 
time dt the point x is displaced into the point % as a result of the 
flow. Then to the first order of approximation the displacement in the 
medium is given by 


(11) x* = x* + u*(x,t)dt. 


3 


It can be shown” that the change 5 of an angle ¢ produced by the 


displacement (11) is given by 


(12) sind 8¢ = [(D, gveveh + D, gush) cost ~- 2D, veveldt, 


to the first approximation, where VY, and Vv, are unit vectors determining 
the angle P and 
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6 
(13) Dig * tag - 3 5B , leap . 3(u, 2 + ap.) Gs ig! 


The tensor D defined by (13) will be called the distortion tensor. 

If ’} = 0 for arbitrary directions V, and Vy, any volume V will be 
displaced into a volume V without change of shape, or as we shall say, 
without distortion. It can easily be shown that the necessary and 
suflicient condition for the displecement (11) to be without distortion 
to the first approximation is given by the vanishing of the distortion 
tensor’. Hence we see from (12) that the tensor D at any point P can be 
considered to give a measure of the rate of distortion, due to the flow, 
of any configuration in the immediate neighborhood of P. 


6. The viscosity tensor. The viscosity tensor V is defined in 
terms of the stress tensor T as follows 


(14) Vig = hg t Poa, 


where p is a scalar function given by p = -7,,/3. The function p is 
called the pressure. 

Py definition a fluid is called ideal if, and only if, the viscosity 
tensor vanishes. For ideal fluids we therefore have "a" Pda and 
hence 7, = -pv,; hence the force T is normal to every surface element. 
In the general case the tangential component to a surface S of the force 
Tis called a shearing force. It follows that for ideal fluids the 
shearing force is zero over arbitrary surface elements and conversely 
it is easily seen that this condition implies the vanishing of the 
viscosity tensor. Hence an ideal fluid can also be defined as one for 
which the shearing force vanishes over any surface element independently 


of its orientation. 


7. Relation between the distortion and viscosity tensors. We now 
make the fundamental assumption that the components of the viscosity 
tensor are linear and homogeneous functions (independent of the time t) 
of the components of the distortion tensor. Thus 


(15) Vag = C22 (x) D,,. 


Since V and D are tensors these relations will be invariant under 
coordinate transformations if we define the quantities Cr} to be the 
components of a tensor. The tensor C so defined, can be taken to be 
symmetric in its two upper and also in its two lower indices due to the 
symmetry property of the tensors D and V. 


8. Homogeneous and isotropic fluids. The equations (15) will be 
specialized by further assumptions. First we assume that the fluid is 
homogeneous. This means that the functions CoB are unchanged by an 
arbitrary translation of the rectangular coordinate axes. Then under the 
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have 
ly —-c) 


Since the c* are arbitrary constants we see from the equality of the 
first and last members of these relations that the values of the C’s are 
independent of the coordinates; conversely if the C’s are constants 
relative to rectangular coordinates the fluid is homogeneous. Hence the 
requirement that the fluid is homogeneous means that the coefficients 
C in (15) are constants in a rectangular coordinate system. 

Next assume that the fluid is isotropic. This means that the values of 
the C’s at any point must remain unchanged under proper orthogonal 
rotations of the coordinates about this point. A tensor possessing 
this property is called an isotropic tensor. It can be shown that the 
most general (symmetric) isotropic tensor C has components given by 


~7J 7) 


(16) 6; 


af 

relative to rectangular coordinates, where A and u are scalars”. When 

the fluid is both homogeneous and isotropic A and “ will be constants. 
Substituting the values of the C’s from (16) into (15) we have 


y = ii 
VB Deo ° 5B . 2uD | 2 ¥ 


When we further substitute for the D,, and the V.2 the values of these 


components given by (13) and (14) we obtain 


(17) 


'aB 


Hence (17) is the most general expression for the quantities 7, 1n a 


homogeneous and isotropic fluid (for which the fundamental reJation (15) 
is assumed). In (17) the scalar 4 is a constant and is referred to 
as the viscosity constant or simply the viscosity of the fluid. Observe 
that if we contract the indices a and / in (17) we obtain the previous 


relation 7,, + 3p = 0 by which the pressure was defined. 


9. The Navier-Stokes equations. We have seen that the equations 
(2) and (3) are equivalent to (10) in which the 7,, are symmetric. 
We now consider the explicit form of (10) when the Tg are given by 


a B,B 


Substituting this expression for the 7, into (10) we obtain 
¢ A 


(18) 
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Making the substitutions (10) and (17) and using (6) this expression 
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and u is taken to be a non-negative constant (the 


where y= u, ~o 
Equations (18) are the equations of motion of a homogeneous, 


isotropic viscous fluid; they are usually referrec to as the Navier- 
Stokes equations. These equations together with the equation of contin- 
uity (5) and an equation of state ($1) furnish us with the theoretic 
means of solving a rather general class of problems of fluid motion. 





viscosity). 













10. Special motions. It is customary in practical applications to 
narrow the class of admissable motions by making certain additional 
simplifying assumptions which can be justified in greater or less degree 
hy the nature of the problem under consideration. Thus in certain 
circumstances we can assume that u = 0 and/or P = const. Also it is 
frequent ly assumed that the velocity vector is the gradient of a scalar 
function in consequence of which we can write ai ee in this case 
it is customary to say that the velocity eduits a potential yw. It is 
usually simpler, moreover, to confine our attention to steady problems, 
i.e., those problems of fluid motion for which the pressure, density 
and velocity are independent of the time. Of the various possibilities 
we note here that the diflerential equations for the steady motion of 
compressible fluids (gases) in the absence of viscosity are the following 


















(19) i? Pe Hat Fs VE F- 










If we suppose, in addition, that 9 = const. and that the velocity 
admits a potential W the equations (19) become 











“ $(ou,u, ) wi ? ws . 


,a 







(20) Rs 










These are the diflerential equations of classical hydrodynamics. But the 
first set of (20) can be integrated and hence can be regarded essentially 
as an equation defining p when the velocity and density are known; in 
this case the problem of determining the fluid motion reduces to the 
solution of the Laplace equation Aj = 0 subject to the given boundary 
conditions. 














II. IDEAL GASES 








ll. The energy equation. The work done per unit time on a moving 
volume V of the fluid by the stresses acting over its surface S is 


given by 





BdS = } 
= J Tapuv dS = J) (ragu®),, dh 








utdV + it dV. 
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for W can he put in the form 


d 
(21) =— | touudV - | p dV + 
dt yy 2 Novo y Moo V 


where ~ is a non-negative function given by 
@ = 2uDig Pap - 


Equation (21) is called the energy equation and the function ¢? is 
referred to as the dissipation function. The first term in the right 
member of (21) gives the rate of increase of the kinetic energy of the 
moving volume V. The second term gives the rate at which work is done 
by the pressure in changing the volume of V. The remaining term 
represents the rate at which work is done on the volume against the 
viscous forces; this work is lost as mechanical energy and must there- 
fore appear in the form of heat developed in the volume V. 


12. The ideal gas. By an ideal gas is meant a gas which satisfies 


the esuation 


(22) 


where T is the absolute temperature, m the molecular weight and RF an 
absolute constant which is the same for all gases. If c, and c, are the 
specific heats of the gas at constant pressure and constant volume 
respectively, then between the quantities R, a, Cp and c, we have the 


relation 
(23) = Jc, - e.), 
where J is the mechanical equivalent of heat. It is assumed that c_ and 


c, are constant for an ideal gas. Defining the constant Y as the ratio 
Cy/ey the above relations (22) and (23) can he written 


p R 
(24) —e Je (y-1)T, -= Je (y-1). 
p v 


m 


13. Thermal conductivity. The total energy of a gas consists of 


its kinetic and internal energies. For an ideal gas it can be shown 
by thermodynamical! considerations that the internal energy per unit 
mass 1s expressed by Jc, T + const. where the constant depends on 
the selection of the zero point of energy. Hence the total energy F per 


unit mass is given by 


(25) E= gu,u, + Jc,T + const. 
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Using (25) the expression (21) for W can be given the form 


4 : { pEdV - J | id d\ dV + { ddV 

, es y a hareeit = = ? P 

_ dt -v" Cv Sy” at yee v 

since the constant appearing in (25) gives no contribution to # in 


consequence of the equation of continuity. 
Now let W' denote the rate at which heat is conducted into the 


moving volume V through its surface S. The value of 4’ is determined 
on the assumption that the rate of conduction of heat into V through an 
element of area dS is given by Jk(dT/dn)dS where k is a constant 
and n denotes distance measured in the direction of the outward normal 
vy to S. The constant k is called the thermal conductivity of the gas. 


Hence 


= Jk J T,v* dS = Jk J Tae dV. 


Now # + W' is equal to the rate of change of the total energy of V 
so that 


d 
(28) — {| pEdV= W+ WW’, 
dt V 


Hence substituting into this relation the expressions for W and W’ given 


by (26) and (27) we find 


dT 
Je | — dV + dV ~ | ddaV - Jk | T dV = 0, 
“v _" dt J pore “y 4y (170 


or, since the volume V is arbitrary, we have 


dT 
(29) Je, p qt © ftigy © Ti gg * O 
This is the equation governing the distribution of temperature in an 
ideal gas. From the equations (22) and (23) we have 


RT 
Je. * Jé.7 «~~ * Jet «=, 
v p m Pp p 
Hence diflerentiating with respect to t and substituting’in the left 
member of (29) we obtain 


dT dp 


= JkT + >, 


(30) J pul 
~~ Za dt e 


when use is made of the equation of continuity. Equation (30) is an 
alternate form of (29). If the viscosity u and thermal conductivity k 
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are assumed to vanish equations (29) and (30) become 
] 


(31) * Se 


14. Entropy. From the discussion in $11 and $13 we see that the 
rate of increase of the heat content of the moving volume is given by 


dS + | pdV = Jk} T,, dV + | pdV, 
V ia V 
Hence JkT _, + D gives the rate of increase of heat per unit volume 
at points in the gas. But this expression is equal to the left member 
of (29) and hence when we divide the left member of (29) by / we obtain 
an expression for the rate of increase of heat per unit mass. This 
latter expression can now be modified as follows. We have 


Jc —+ 


"at 


ta [dT read 
7 “v lat r= 0 dt | 


dlogT dl 
JeT js tes? o tpoes BO, 
7 ee dt 


d Pp ) d logp 
Je T |}—1 —|- (y-1) ——], 
"° Fe oe ; dt | 


d p } 
YS ee 5 le 5 ‘ 
”. 2 p 


O,o O,o 


dT p 
p 


dT | 
u =Jc |— + (y-1)Tu 
° dt il 





4 


We now define a function? by the equation 
p 
(32) T=Je loe| 2+ const., 
v | p” 


where the constant is absolute, i.e. independent of position and time. 
Hence 7 dz/dt is the rate of increase of heat per unit mass. The 
function 7 is called the entropy per unit mass or simply the entropy of 
the gas and for any particular gas it is given, by definition, to 
within an additive constant. Solving (32) we can also write 

elev 


(33) p= N(o)p”, N(x) =C 


where C is an absolute constant for any gas. 
If the thermal conductivity k and viscosity / are both equa! to 
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zero it follows from (29) that do/dt = 0. Hence in an ideal gas with 
conductivity and viscosity zero the entropy is constant following each 
particle. In particular if the motion is steady the entropy is constant 
along each stream line; this implies that the function N (a) in the 
relation (32) is constant along each stream line although this func- 
tion can vary from stream line to stream line. When the motion is such 















bs that NV is constant throughout the gas we say that the gas satisfies the 
adiabatic condition or that the flow is adiabatic. The assumption that 
the flow is adiabatic results in considerable simplification of the 
mathematical analysis but is, strictly speaking, an untenable hypothesis 

- except in the case of certain especially selected problems. 

le 

ber 15. An equivalent of the entropy condition. Assuming k = u = 0 so 

ain that do/dt = 0 let us differentiate (33) totally with respect to t; on 

lis making use of the equation of continuity (5) and the equations of 





motion (18) with 4 = 0 the resulting equation can be written 











3 
(34) - - Pus —— — PU, BU, Ug * Tite, 








Conversely from (34), the equation of continuity, and the equations of 
motion with = 0 we can show that doz/dt = 0 in the relation (33). 
Hence in an ideal gas with thermal conductivity and viscosity equal to 
zero the condition do/dt = 0 is fully expressed by (34). In case the 
motion is steady the condition do/dt = 0 is given by 









(35) Pu, g4.4g - Ypu,, * 0. 









16. The differential equations for ideal gases. The complete set of 
diflerential equations by which the motion of an ideal gas is determined 
is given by the equation of continuity (5), the equations of motion (18) 
and the equation (29) for the distribution of temperature in the gas; 
in addition there is an equation of state (22) which was originally 
introduced as the defining condition of the ideal gas. It is of course 
possible to eliminate the temperature T between (22) and (29). When 
this is done we have in place of (29) a differential equation in the 
quantities p, 9 and u,; the relation (22) then appears as an extraneous 
equation defining the temperature 7. In particular the steady motion of 
an ideal gas without viscosity and thermal conductivity is governed by 
the differential equations, 





















(36a) 





(pu,) , = 0, 











(36b) Ra * Py Sym * 0, 





(36c) Pu, guaug - Ypu,, = 9. 
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As we saw in $15 the last of these equations is equivalent to the 
entropy condition do/dt = 0 and hence from $14 it is equivalent to (29) 
with k = w= 0 or to either of the equations (31). 


17. The Bernoulli equation. Multiplying (36b) by u, and summing 
on the repeated index a we obtain 


dp 


d 
—+ $p — (u 


dt dt atta) = 


Dividing this equation by ( and making use of (33) in which o = const. 
along the stream lines we find that 


(37) 


+ su u, * COnSt., 


Y p 
y-lp 


where the constant depends on the stream line. Equation (37) is the 
Bernoulli equation. It holds for the steady motion of an ideal gas for 
which the viscosity and thermal conductivity are zero. 

When we assume that the velocity admits a potential ($10) equations 
(36b) can be written 

Pia. S$(uu). = 0. 
CGC 8 ,a 

If, in addition, the flow is adiabatic ($14) the pressure p in these 
equations can be eliminated by the substitution (33) and the resulting 
equations can be integrated. This leads to an equation of the form (37) 
in which the constant in the right member is independent of the stream 
line. To distinguish between these two cases the equation (37) with 
constant depending on the stream line is sometimes referred to as the 
weak Bernoulli equation and when this constant is absolute, i.e., 
independent of the stream line, the equation is called the strong 
Bernoulli equation. 


III SHOCK WAVES IN IDEAL GASES 


18. Singular surfaces. Consider a moving surface 2(t) such that 
each point of >(t) for t fixed, lies in a regular surface element! and 
let f(x,t) be a function which is discontinuous in the variables x° at 
all points of S(t). For points x* not on > at any time t we assume 
that f(x,t) is continuously diflerentiable in all variables; it is 
assumed also that f(x,t) and its derivatives exist and are continuous 
on each side of 2. The surface >(t) will be said to be singular relative 
to the function f(x,t). 

Now let V be any moving volume in the fluid which is divided by the 
above surface 2 into two volumes V, and V,. Let S be the sur face 
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of V and denote by S, and S, the portions of S which form parts of 
the boundaries of V, and V, respectively; the remaining part of the 
boundary of V, and V, will then be furnished by the surface 2. As 
in the foregoing discussion we suppose that the variation in the volume 
V is produced by the moving particles of the fluid so that the normal 
component of the velocity of V at points of its surface S is given by 
u, = u,V¥* where v is the outward unit normal vector to S. Let G 
denote the normal velocity of =(t) along the outward normal to 2 consi- 
dered as part of the boundary of V,; then -G will be the corresponding 
normal velocity when 2 is taken as a part of the boundary of V,. Now 


8) — | f(x,t)dV = — J f(x, t)dV — { f(x,t)dV 
) ;* in ,t ij j ,t ® 


dt ; V, 


But each term in the right member of this equation can be evaluated in 
accordance with the formla (4). This gives 


d of ) 
— | f(xz,t)dV = — dV+j} fu,dS + j_f GdS, 
V y ot S + 
1 ' 1 

no 1 a 

f(x,t)dV = —dV+!} fu,dS - j_ f,GdS, 

V, J, 3 ‘Sp 5/2 
2 

where f, 1s the value of f on the side of 2 bordering V, and f, the 
value of this function on the side of 2 bordering V,; the designation 
2 as the area of integration in the last integrals in these equations, 
as well as in corresponding integrals in following equations, naturally 
refers to the part of the surface 2 bounding the volumes V, and V,. 


Hence substituting into (38) the resulting equation can be written 


of 


ot 


d 
(39) aati y = 
— J fla, td i 


dV + J funds + Lf, - f,)GdS. 


19. Shock waves. We now assume that the surface 2(t) is sin- 
gular relative to the pressure p. The quantities P and u, will be 
assumed to be continuous and differentiable for points x not on 2(t); 
also these quantities and their derivatives are assumed continuous on 
each side of 2(t). The moving volume V, which is divided into the 
volumes V, and V, of $18, must satisfy all integral conditions previ- 
ously imposed and which have resulted in the differential equations 
for the motion of the gas under the hypothesis of complete continuity. 
These integral conditions are the following: first, the law of conser- 
vation of mass, second, the assumptions (a) and (8) of $1 and third, 
the assumption contained in the equation (28) which led to the differ- 
ential equation for the distribution of temperature in an ideal gas. By 
allowing the volume V toshrink to zero we thus obtain certain relations, 
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called shock conditions, between the values of the quantities p, ( and 
u, on the two sides of the surface 2(t). A surface >(t), singular 
relative to the pressure p, over which these shock conditions are 
satisfied is referred to as a shock wave. In the following sections we 
derive the shock conditions for ideal gases for which the viscosity and 


thermal conductivity are zero. 


20. First shock condition. For convenience of terminology we shal] 
refer to the side of > bordering V,; as the side 1 and to the side 
bordering V,as the side 2 of 2. Then from (1) and (39) with f = p(x,t) 


we have 


ra dv + { dS d lL. ( )Gds 0 
ann + 5 + o> i. - dS = 
ES . J.P Un Js \P1 — Pe 


where (, and P, denote the value of the density along the sides 1 and 2 
of >. Now let V approach zero at a fixed time t in such a way that in 
the limit it passes into a finite part 29 of the surface >. Then the 
volume integral in the above equation is of higher order than the 
surface integrals and can be neglected in comparison with the latter. 


Also® 


J pu,dS~ - f. pu, ds, 
Ss; 4 2 | 


f pu,dS > f Py Uz,dS, 
Ss." . 


where u,, and uz, denote the normal components of the fluid velocities on 
the sides 1 and 2 of = along the normal direction from the side | to 
the side 2. Hence we obtain 


Js (0, Ur, — GP, )dS “Js (P2 Ury - Gp, )dS = 0. 
0 


Since this condition is independent of the extent of the surface of 
integration 2), we have 


(40) P, (u,, - @) = p, (u,, - G) 
over >(t). This is the first of the shock conditions. It originates 
entirely from the principle of the conservation of mass. 


21. Second shock condition. The assumption (a) of $1 has its 
expression in the equation (2) in which we suppose the applied force 
X = 0. Assuming the viscosity “ = 0 so that T,= -pv, the condition 
(2) thus becomes 


(41) TS, puadV + J. py, dS = 0. 





HYDRODYNAMICAL EQUATIONS FOR GASES 183 


We now allow the volume V to approach zero as stated in $20; it is then 


easily seen that 


0 
where the Vv, in the limit integral are the components of the unit normal 
directed from the side 1 to the side 2 of >. Also applying (39) with 
f =u, the first term in (41) becomes 
‘a 
ee } { . { » . 
oe + J, Puan ds + 1g, Pa + I Oy Uy, P,u,, )GdS, 


and this expression, on passage to the limit, yields 


(43) Js P,u,,(G-u,, )dS - Jo Pz Uzq(G- up, dS. 
“0 ~O 

Hence substituting (42) and (43) for the two terms comprising the left 
member of (41) and making use of the first shock condition (40) we are 
led to the relations 
(44) [p] yu =p, (G-u,,) Cu), 
where we have put [p] = P>-P, and Cu] = U,,-U,,- The brackets ( ] 
are here used to denote the change in a quantity in the passage over the 
surface >. This conventionwill be adhered to in the following discussion. 

Turning now to the assumption (8) of $1 which is expressed by the 
equation (3) it can be shown that no new relations are obtained beyond 
those given by (40) and (44). Hence (44) constitutes the second of the 
shock conditions according to the enumeration of $19. We emphasize that 
this shock condition holds for fluids or gases for which the viscosity 
is equal to zero. 

The product p,(G-u,,) cannot vanish at any point of the surface 
2. For if this quantity would be equal to zero it would follow from (44) 
that [p] = 0 thereby contradicting the condition for the surface = to 
be a shock wave ($19). The fact that the product p,(G-u,,) does not 
vanish, and hence that 0,(G-u,,) does not vanish by the first shock 
condition, is used in the following derivations. 

Multiplying (44) by v, and summing on the repeated index a we can 
deduce the relations 


(45) — [p] = o,(G-u,,)lu,) 


by which [ug] and [p] are expressed in terms of [u,]. Conversely (45) 
implies (44). Hence (45) is equivalent to the shock condition (44). It 
follows from the second equation (45) and the condition [p] # 0 that 
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{u,] can not vanish at points of a shock wave >; using this result the 
first equation (45) implies that not all of the quantities [u,] can he 
equal to zero at any point of 2 if this surface is a shock wave. 


Again, using (40) and (45) we can deduce the relations 


(46) ete shes 
js 





) ; 2 )] 
(47) (p] Py (G~uy,) lol 


P. 
from which the conditions (40) and (45) can be recovered. Hence the 
conditions (46) and (47), expressing [u,] and [p] in terms of [p], are 
equivalent to the first and second shock conditions. 

Finally observe that if we multiply (44) by A*, where 
vector tangent to >, we obtain [u,r*] = (0. In other words the tangential 
components of fluid velocity are continuous across the surface 2. 


» 
\ 


is a unit 


22. The Rankine - Hugoniot relation. We now impose the ful] 
restrictions mentioned in $19, namely that we are dealing with an ideal 
gas with viscosity and thermal conductivity equa] to zero. Then the 

° ’ Nd e S ° e 
quantity # of $13 vanishes and 7, = -pd_,; hence substituting for 4 
ap ap 


in (28) this equation hecomes 


d : 
( 48) —— pEdV + pu, v ‘aS = (). 
dt “V S 


But from (39) with f = OE the first term of this equation can be written 


(OE )dV + OEu,dS + | 1(E) 


S > 


17 (0E) ,}GdS. 


Now substitute this expression for the first term in (48) and then allow 
V to approach zero as in the preceeding discussion; then making use of 
the first shock condition (40) the resulting equation leads to the 


condition 
[ Lae rey 
(49) pu, | p,(G-u,,)(E]. 


If we put i = Jc, T it follows immediately from (22) and (23) that 
Je,T = .-p/o. Making this substitution for Jc, T in the expression (25 
for E we see that the relation (49) can he written 


| . r >] 7 
(50) [pu,] + p,(G-u,,) j 2 | -p,(G-u,,)li+ su, ug] 
LP J 


Now using (40), the second relation (45), and (47) it can readily be 
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shown that 






[pu,] = {o,(G-u,,)u,, + py}lu,), 


al 


Then when we substitute these expressions for [pu,] and [p/e] into (50) 


we obtain 
[i+ su,u,] = Glu, ] ; 


The relation (51) is one of the forms of the third shock condition. 
From the relations of $12 and the conditions (46) and (47) it follows 


that 
y k 
y-1Le]} ’ 


Hugi] = Gla) - | 





P, [u,] 


. p,(G-u,,) 


(G-u,,)lu 





n 







(51) 






[t] = 








P,*P2 
2p; P2 





|). 






Making these substitutions in (51) the resulting expression can be 


written 






P2 _ (y+1l)p, - (y-1)p, 
Py : (y+l)p, - (y¥-DV)p, 





(52) 












The form of the third shock condition given by (52) is usually called 
the Rankine -Hugoniot relation. It holds for an ideal gas in which the 
viscosity and thermal conductivity are zero. It is an interesting fact 
that the form (52) of the third shock condition does not involve the 
wave velocity G even though it is not confined to steady shock waves. 
The condition (52) can be written in the following equivalent form 















(53) Po 7 (y+1)p, i (y-1)p, 
Pp, (Y¥-l)p, + (y+ 1p, 








Putting x = ~,/o, and y = p,/p, the above relation becomes 





(y+l)y + (y-1))) (y#1) Ay 
~ (yeldy e (¥4+1)) yel— (y-1)?y + (y? -1) | 









Since Y> 1, x will have its greatest value when y= © and this maximum 
value of x is (y+1)/(y-1). The least value of x occurs when y = 0 and 
ls given by (y-1)/(y+1). Hence 






THOMAS 


For air Y = 1.4 approximately; for this value of y the above inequality 
becomes 1/6 $ x § 6. 


23. The combined shock conditions. The shock conditions which we 
have deduced are expressed completely by the relations (46), (47), and 


(52). It is a simple algebraic consequence that these relations can be 
put in the following equivalent form 


y+ ] 


9 f a © 
AP loot hei | 


2yp, - lyn Dy AG=— ty 


by which the values of the quantities u,, p and on side 2 of the 


surface > are given directly in terms of the u,, p and ( on side | of 
the unit normal Vv to 2 (directed from side 1 to side 2) and the velocity 
of propagat ion G of the wave surface. In the case of a stationary wave 


front for which G re relations (54) hecome 


w «¢ + ee 
From the remark in $21 it follows that the denominator in (54a 


CAL) 
cannot vanish at points of a shock wave. Also the numerator in (954b) 
cannot vanish since otherwise we would have |p| = 0; hence the denomina- 


tor in (54c) must be diflerent from zero as its vanishing would imply 


i 


[ep] = w (see $19.) 
24. Change in entropy and density in passage through a shock wave. 
The thermodynamic condition da/dt 2 0 for continuous flow without ther- 


mal conductivity leads to the requirement that the entropy per unit 
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mass 7 associated with a moving particle does not decrease in the 
passage of the particle through a shock wave (considered as representa- 
tive of the rapid variations in a flow which is actually continuous). 
Thus if the particle moves from side 1 to the side 2 of a shock wave 2 
we must have 7, 2 9, and this implies p,/Pz 2 p,/P?} on account of the 
expression (32) for 7; hence 


Pe 
me) 60 a 
Py Py 


Putting x = 0,/0, and using the Rankine —-Hugoniot relation (52) for 
p,/p, the inequality (56) becomes 


1) (y+ 1l)x - (y-1) — 


(y+1) = (y-B)s * ” 





Let us now define a function y by the equation 


(y+1l)x - (y-1) " 
(y+1) - (yv-l)x 





(58) y = 


Then y = 0 for x = 1; also both dy/dx and d*y/dx* vanish for x = 1 while 


d3y 
—3 = 4(y-1)(2y-1)3-y). 

dx 
Taking Y = 1.4 we see that d°y/dx*® > 0 for x = 1 so that the graph of 
(58) has a point of inflection at x = 1. The following table gives approx- 
imate values of y for the indicated values of x. The graph of this 
function is shown below. Condition (57) is satisfied only for values 


y 
-1/6 
-.08 
0 
11 
1.01 
4.54 
19.48 








Nn 


























OY em] WiDr O]& 


of x for which y 2 0. Hence we must have ~,/0, 2 1. If Po = P, the 
relation (47) gives [p] = 0 contradicting the assumption that > 
2 is a shock wave. It follows that P> > P,;. But x > 1 means that y > 0 
and hence the inequaltiy sign holds in (57). This implies 7, > 7,. 
Hence when a fluid particle traverses a shock surface the density and. 
entropy per unit mass associated with the particle are increased. 

If a gas increases in density when it passes through a shock wave, 
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the wave is called a compression wave. Hence shock waves, as here 
defined and treated, are compression waves in view of the above ther- 


modynamic condition. 


NOTES 


1. O. D. Kellogg, Foundations of Potential Theory, J. Springer, Berlin, 1929. 
See p. 105 for definition of a regular surface element and p. 112 for the 
definition of a reqular surface. 


2. The usual proof of this relation involves the assumption of the equality of 
action and reaction in consequence of which the stresses T; and T2 on the two 
sides of a surface are equal in magnitude but opposite in direction. However 
this can be shown by use of the condition (2). For let T; and T> be the 
stresses on the two sides of a regular surface element > (see footnote 1, loc. 
cit. p. 105). Let V be a small sphere which is divided by the surface element 
into volumes V; and V5. Then applying (2) we have 


d 
(a) ey Puad¥ = J, TadS, 


d f f 
—— | Pagd¥ * | TadS * J.%, 05, 
dt A S; > 


d ( f 
(c) +f punav=f{ T,aS + { %odS, 
Cc de J V, Ua ‘5, a y 2a 


where S, and S, are the portions of S bounding V, and V, respectively and it is 
assumed that the stresses are continuous over the surface S of the sphere and 
on each side of the element 5. But adding corresponding members of (b) and (c) 
and using (a) we obtain 


Ene Deas = 0. 

Js Na Js Tea 

Since the area of integration 2 is arbitrary it follows that T,q + To, = 0. 
The remainder of the derivation of the relation (7) can be made by considering 
the stresses over the surface of a tetrahedron three of the faces of which are 
perpendicular to the coordinate axes. See P. Appell, Mecanique Rationnelle, 


vol. 3, 1928, p. 134. 


3. Cf. A. E. H. Love, The Mathematical Theory of Elasticity, Cambridge, 4th 
ed., 1927, p. 62. Also S. Timoshenko, Theory of Elasticity, McGraw-Hill, 1934, 


b. 392. 


4The vanishing of the distortion tensor is evidently sufficient for 5¢ = 0 in 
(12). To show that this condition is also necessary we have 


(a) (Dagurry + Dapreve) cos gd - Da pry = 0 

for arbitrary directions V,; and Vy. Taking V,; perpendicular to Vy so that 
cos ? =0 equation (a) gives Dp = 0. When Y has the direction of the 
x? axis and Vp the direction of the x” axis, where 747, the latter equation 


gives D,,,=0 (a4). Nowchoosev, sothat this vector does not make a right angle 


with any of the coordinate axes. Then taking v, in the direction of the positive 
x’ axis we have cos ¢ = vive = v} 7 0. Under these conditions (a) becomes 


a 


a 


In a similar manner by taking Vv, in the directions of the positive x and x 
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axes, but leaving the above vector , unchanged, we obtain 


(c) DapvaUs = Doo, 


(d) Dae) = D33- 
Adding (b), (c), and (d) we find that 
(e) Dapsve = 0, 


since D,, = 9 from the equation (13) defining the deformation tensor. Now since 
the expression in (e) is homogeneous in the quantities this equation is seen 
to hold for arbitrary values of the 3, i.e., independently of the original 
condition that V, is a unit vector. Hence (e) implies Dag = 0 or, in other 
words, the vanishing of the distortion tensor. 


5. H. Jeffreys, Cartesian Tensors, Cambridge, 1931, p. 68. 


5. The first of these limits can be inferred as follows. Under the conditions 
of diferentiability and continuity assumed we have 


j (Pua). dV = j PundS + | Py un dS. 

V, S; 2 
Then allowing V to shrink into the area 26 the left member of this equation 
approaches zero and hence the limit in question is obtained. The other limit 
can be inferred in a similar way as well as the corresponding limits in the 
following sections. 


T. Y. Thomas 


Indiana University 
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Abstract: The relation of isomorphism is investigated at some length 
and after consideration of the usual uses of the term and its associated 
notion generalization is developed. An abstract isomorphism is defined 
so as to contain the usual sense on specification. Some of the aspects 
of this consideration are new therefore some positive contribution is 







made. 
As an aid in thought suppose I arrange a collection of piles of 


pebbles on the table. Actually you wil] have to imagine my performance 
of the act since I have no intention of undertaking the arduous labor 
of creating such an object assembly. Instead I shall look at the whole 
deal abstractly and place a collection of dots on the blackboard. 













] 2 3 4 bs) wre 13 . 









The dots are not pebbles, they merely represent pebbles. In a 
different context I might represent an element of a group by a matrix. 

Further I will perform an abstraction and represent each represen- 
tation of a pile of pebbles by number symbols. In the case of the pile 
of thirteen pebbles I have used thirteen dots and represented these 
symbols by the configuration of symbols 13. As I write the symbols I 
simultaneously utter an auditory representation which represents the 
language representation and does not need translation to you. At least 
we agree there is a lot of representation going on around here. 

The integers are no mere trite representation. They were obtained by 
unconscious operations of intellection, representation and abstraction, 
and they have the good fortune to be universals which are indispensable 
to thought on any level. Each integer represents a collection of classes 
of objects. 

To further consider the notion of representation we could consider 
all the things a man represents. That, however, would lead us astray 
from the subject of this paper. Such modes of analysis are very useful 
and many a false argument can be laid bare with a minimum of effort. For 
example, Prof. Max Black in Principles of Really Sound Thinking which 
appeared in the March 1948 Scientific Monthly gave as maxim number three 
a all depends on who says it.’ In support of this maxim he says: 

Men are easier to classify than arguments—attend to the man, not the 
argument.’’ Such a slight argument hardly needs refutation since one 
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only needs to think of the order of complexity of the objects under 
contemplation. A man, any man, is more complex than any argument there- 
fore Prof. Black’s argument is false. 

To return, the representation on the blackboard is one-one. A set 
of dots corresponds to an integer. The collection of sets of dots 
corresponds to the collection of integers, one-one. Also the hypo- 
thetical collection of sets of pebbles stands in biunique, one-one, 
correspondence with the collection of sets of dots. 

If we add the set of two dots to the set of three dots we get a set 
which is equivalent to the set of five dots. In a corresponding way when 
we add the integers two and three we get five. Furthermore we have 
assembled the parts of the notion of isomorphism. The collection of sets 
of dots is isomorphic to the collection of integers; that is, there is 
a one-one correspondence between them and under this correspondence the 
operation of addition is preserved. This is the basic meaning of the 
term isomorphism. 

Thinking in a different mode and proceeding from automorphism 
MacDuffee defines isomorphism and then states ‘‘Two isomorphic math- 
ematical systems are abstractly identica!.’’ Such a unique statement 
deserves repetition here because each word contains so much meaning. 

Isomorphism is a true equivalence relation. While demonstrating this 
instead of the usual symbol = we will use a more general R. If we let 
a represent a typical integer, 6 the corresponding set of dots and c the 
corresponding pile of pebbles we can give the standard demonstration, 

Re flexive: aRa for all a 
Symmetric: aRkb implies }bRa, 
Transitive: afb and bRc imply akc 

When isomorphism is applied to different mathematical systems special 
conditions have to be met. For two domains to be isomorphic not only 
addition but multiplication too must be preserved in the biunique 
correspondence. The same holds true in the case of fields. Our defini- 
tion is closest to the case of groups where only one operation, the 
group operation is involved. This group operation must of course be 
preserved and this One condition is the sole difference between an 
ordinary one-one correspondence and an isomorphism. 

What I propose to do is to abstract a generalization from the iso- 
morphism as defined. Instead of saying addition is preserved I wil! say 
the isomorphic operation remains invariant. 

I define abstract isomorphism to be a relation containing a biunique 
correspondence under which a set of operations is invariant. If the set 
of operations is a null set the abstract isomorphism is merely a one-one 
correspondence. The question arises as to whether such a case exists. 
I am still investigating the structure of operations so that a statement 

can be made covering all operations. Even specification by which the 
usual senses of isomorphism are shown to be contained in this definition 
is itself an operation so I’]] have to restrict my remarks to a simple 
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demonstration of inclusion. 

If | specify that the set of operations consists of addition and 
mltiplication I readily obtain the sense of isomorphism as it exists 
between two domains or between two fields. If I specify that successive 
substitution is the set of operations the term abstract isomorphism 
takes the sense of isomorphism between substitution groups. Since 
under this new definition no operation is required it can be said that 
an abstract isomorphism exists between matter waves and electromagnetic 


waves. 
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ANNOUNCEMENT OF TWO SYMPOSIA ON NUMERICAL ANALYSIS 


The National Bureau of Standards is planning two symposia on the 
effective utilization of automatic digital computing machinery to be 
held in June 1949 at the Bureau’s Institute for Numerical Analysis in 
Los Angeles, California. 

Symposium I: Construction and Applications of Conformal Maps. Appli- 
cations of conformal maps in such fields as aerodynamics and electronics 
will be emphasized, with special attention devoted to the current 
needs of research workers. ° 

Symposium II: Probability methods in Numerical Analysis. The second 
symposium will cover probability methods in numerical analysis, with 
emphasis on applications of the “Monte Carlo method” that has proven 
so useful in the solution of many problems in mathematical physics. 


Those who would be interested in attending either of these symposia 


may obtain further information from Dr. J. H. Curtiss, Institute for 
Numerical Analysis, Los Angeles 24, California. 








Announcing two more mathematics texts... — 


ANALYTIC GEOMETRY 


By ALFRED L. NELSON, KARL W. FOLLEY, and WILLIAM M. BORGMAN of 
Wayne University 


REPARED for use in a college course in analytic geometry, this text is planned a 
preparation for the calculus rather than as a study of geometry. In order that it may bh 
of maximum value to the future student of the calculus, the basic sciences, and engineering 


considerable attention is given to two important problems of analytic geometry. They ar 


1) given the equation of a locus, to draw the curve, or describe it geometrically; (b 


iven the geometric description of a locus, to find its equation 215 pages 


INTRODUCTION TO ANALYTIC GEOMETRY 
AND THE CALCULUS 


DADOCURIAN, Trinity C i(lLonn 





By EARLE B. MILLER, Illinois College 


. LEAR, carefully organized exposition written primarily for students who have had 

only one year of algebra in high school. Since its publication a short time ago, we have 
received enthusiastic comments from colleges and universities throughout the country 
“Professor Miller is on the right track. There is @ great need for an intermediate cours 


in algebra.”—Professor Daniel W. Snader, University of Illinois. 361 pages $2.5( 


PREPARATORY BUSINESS MATHEMATICS 


By LLOYD L. SMAIL, Lehigh University 


HIS BOOK was written specifically to meet the need for a text which will give the 

business student adequate preparation for subsequent courses in mathematics. Of this 
text Professor D. H. Lehmer of the University of California has written, “I find the topics 
treated are painstakingly well done, as is the case in all of Smail’s texts. The figures and 
tables are very good, too. It should make a good second semester course following algebra 
and preceding mathematics of Finance for Business Administration students.” 


244 pages $2.75 
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EQUIPARTITION OF CONVEX SETS 


R. C. and Ellen F. Buck 


It is well known that a line can be drawn thru any interior 


perigee ' 
f » convex set dividing the set into two pieces of equal a 


What is Mathematics, Courant and Robbins (p. 317 et 


f the same type are proved. In particular, 1t 
ee eer ‘ 
region’ can always be divided into four equa 


uts, and that furthermore they can be chose 


| 
meceerrnet 


into eve 


“A set is convex if the line segment joining any two of its point 


in the set.” 
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Move P along the boundary until it reaches the position now occupied 
by P’. Since the line PP’ is the only line thru P’ bisecting the region, 
P’ will have moved along to take the former position of P. Again by 
uniqueness, Q and Q’, and R and R’ will have exchanged positions. Let 6 
be the angle RNR’ measured counterclockwise; @ is a continuous function 
of P. The sum of the old and new values of @ is 27; if one of the values 
is greater than P, the other is less. Thus, for some position of P 
between these extremes, 9 will be exactly 7; the ‘ine RNR’ is then 
straight and the region has been divided as desired. 

Where there is general intersection of the three cuts, instead of 
concurrency, the system has another degree of freedom. It might be 
expected that the lines could then be so chosen that all of the resu!t- 
ing regions have equal area. This is not true if the region is convex. 

Suppose we have a convex region of area seven, divided by the lines 
EF, GH, and DI into seven regions each of area one (Figure 2). 


Starting from this assumption, we shall arrive at a contradiction. Draw 
the lines DE, AI and FG. The area of the shaded quadrilateral CBDE is 
less than that of the larger curved part containing it, which is of 
area one. The same is true of the other two shaded quadrilaterals. 
Instead of the convex set we shall consider the triangle PQR and the 
seven parts into which it is cut by the dividing lines. The triangle 
ABC has area one; we can take the shaded quadrilaterals of figure 2 
to have area one, and the remaining quadrilaterals BGPD, CEQH, and 
AIRF to have area not less than one. For, otherwise we move lines DE, 
GF and IH outward by parallel displacements until the areas of the 
shaded regions are exactly one each, and this will not reduce the areas 


of the other three quadrilaterals. ” 
We now observe that of al] the permissable positions of the line FG 
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such that the area of FGBA is one, BGPD will have the greatest area 
when BG is greatest, and hence when FG passes thru A (Figure 3.) We 
shall prove that even in this position, BGPD has area less than one, or 
if it is one, AIRF has area zero, contradicting our assumptions. 


Suppose that ED is parallel to or converges toward CB as in Figure 3. 
Thru D draw E’ DP’ parallel to CBG. Triangles ABC and ABG have area 
one. The quadrilateral DECB has area one, and hence the area of BDE'C 
does not exceed one, so that the area of ADE’ does not exceed two. By 


similar triangles 


area ADE’ AD* _ area AP'D < 


area ABC AB? area AGB 


so that area APD $ area AP'D $ 2, and area BGPD £ 1. The equality 
sign arises only when ED is parallel] to CB; otherwise BGPD has area 
less than one, contradicting our assumption. In the extreme case, FG 
must pass thru A and the area of AIRF is zero, instead of being not 
less than one. (If ED converged in the opposite direction, the argument 
would be the same, using the quadrilateral CHQE, Figure 2). 

This completes the proof: no convex set admits an equipartition 
by three lines. Simple examples show that some non-convex regions can 
be so divided; however, the converse of the theorem is not true since 
examples can also be constructed which are not convex but which share 


the property of admitting no such equipartition. 
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We have actually proved more than the above. If in Figure 2, the 
curved pieces of the convex set labled BCED, BDG, GFAB and the central] 
piece ABC have area one, then the piece FAI has area zero. As an 
immediate corollary of this, if a convex set is divided into seven 
pieces, of which six are equal, the seventh (unequal) piece must be 
the central piece, ABC Figure 2. If the equal pieces have area one, 
there is a least upper bound for the area of the seventh piece; we 
conjecture that the seventh piece can have area 1/8 at most. This 
bound is obtained by dividing a triangle with cuts parallel to the 


sides. 
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A NEW LEAST-SQUARES LINE* 
G. A. Baker 


The usual method of fitting a straight line to a set of data as given 
by Snedecor (1) say, is to assume that the independent variable, x, is 
free of error and that the variance of the dependent variable y is 
uniform for all x’s. In considering the relationship between white and 
yolk weights of birds’ eggs, Asmundson et al (2) found that the usual 
assumptions were not sufficiently close to reality. A study of the 
actual data suggested that a close approach to reality would be obtained 
by considering a least-square line with the distances from the observed 
points to the fitted line in some direction weighted inversely as the 
squares of the distances of the points of intersection of the lines 
through the points with this direction with the fitted line from the 
y-intercept of the fitted line. The purpose of this note is to discuss 
fitting a line under these conditions. 

Suppose that we assume a,configuration of points (x, ¥;) in the 
first quadrant and that the least x. > 0. We shal] assume that x, and 
y; are approximately connected by the relationship 


(1) y = ax + B, a # 0. 


We shall consider the distances, d,,, from the points (x., y;?) to the 

line (1) taken in some direction m not necessarily perpendicular to (1). 

We shall consider that the weight to be attached to d,,; 1s inversely 

proportional to the square of the distance, d,s: of the point of inter- 

section of the line through (x;, y,;) with the direction m from the 

y-intercept of (1), which is (0,8). 
We compute 


2 2 
ae 2\2 


(a —m 


2 


-* (a = a)* 


Thus, we wish to determine m, a, and & so that 
2 
L+m (y; - ox; 





2 
(4) f (nm, a, p< De —- + ae (y. —- MX- — PP 
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Part of a paper presented under the test le of, “SA Mathematical Model 
of the Relation between White and Yolk Weights of Birds Eggs, at the 
twenty-eighth meeting of the Institute of Mathematical Statistics, June 


17—19, 1947, San Diego, California. 
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shall be a minimum with variations in m, a, and §. Necessary conditions 


for such a minimum are 


“2 
(CQ + mn )(y, - ax, - §)'x, i (y, - ax, 








(5) 


t 
L (y, - mx, - B)° (y, - mx, 


r e 2) 2 
(y, -ax. - Bix, (y, - ax. - B)*a 
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L(y, - mx, - B) (y, - mx, - B)*(1 + a?) 
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Now, following Deming (3) and finding the least-square line with 
unequal weights (not specified here but based on empirical determina- 
tion) and both x and y subject to error we obtain 


(8) y = 0.762 + 1.715x, . 

We use values from (8) as first approximations to a, 4, and m. That is, 

replace a by 1.715 + a; 6 by 0.762 + 6; and, m by -0.583 + c since 

m= — — approximately where a, 6, and c are now regarded as small. 
Expanding (5), (6), and (7) after replacement in powers of a, b, and 

c and neglecting higher terms we obtain 

(9) a = — 0.155, 6 = 0.206, c = ~ 0.006 

and hence the equation 

(10) y = 0.968 + 1.560x 


as an estimate of the relationship between the two variables white and 
yolk weights of birds’ eggs. The angle between the direction m and this 
‘‘best’’ fitted straight line is 87° 50’ instead of 90° as assumed in 
deriving (8). The fit of the line (10) to the observed data appears to 
be good. 

The conventional .least-square lines measure deviations from the 
fitted line parallel] to the x-axis, parallel to the y-axis, or perpen- 
dicular to the fitted line. The present method generalizes the direction 
along which deviations from the fitted line are measured. A study of the 
data available indicated that the weights chosen which were approx!- 
mately equivalent to the square of the egg size were sufficiently close 
to reality to give an acceptable final result. 
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TEACHING OF MATHEMATICS 
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Joseph Seidlin and C. N. Shuster 


This department is devoted to the teaching of mathematics. Thus articles on 
methodology, exposition, curriculum, tests and measurements, and any other 
topic related to teaching, are invited. Papers on any subject in which you, as 
a teacher, are interested, or questions which you would like others to discuss, 
should be sent to Joseph Seidlin, Alfred University, Alfred, New York. 


ANALYTIC GEOMETRY—THE FRAMEWORK OF MATHEMATICS 
Charles K. Robbins 


The use of graphs dates back to the astronomer Hipparchus and 
probably to the ancient Egyptians. In modern times graphs have been used 
extensively to present statistical and scientific data. As an example 
consider the following table which displays simultaneous readings (to 
the nearest degree) on Fahrenheit (denoted by y) and Centigrade (denoted 
by x) thermometers. 


x -20 -10 0 10 20 30 40 
y ~4 14 32 50 68 86 104 


To present these data graphically, the centigrade readings are 
measured on a horizontal line as distances from some point O, positive 
to the right of O and negative to the left (Figure 1). At x = 10 we draw 


/' 


/ 





/ 
































Figure 1 Figure 2 


a line 50 units long vertically upward, at x = 20 a line 68 units long 
and so on. At x = -20 the line is drawn 4 units downward. We thus locate 
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7 points. If we were to take further readings for x between 10 and 20, 
say, then the readings for y would be between 50 and 68. Therefore 
continuous curve drawn through the 7 points represents, at least approx- 
imately, the relation between x and y. 

Often valuable information is made obvious by a graph. In this case 
the graph appears to be a straight line. If it is a straight line we 
say that the relation between x and y is linear. Also from the graph we 
can read (appreximately) a value of y corresponding to any value of x, 
Thus if x = 5, y = 41; if y = 0, x = -18; and so on. 

The study of graphs leads us to generalize the process of locating 
points in a plane with reference to any two mutually perpendicular 
lines. If you stand in a fence corner and say: ‘‘I am 5 yards east of 
the North-South fence and 3 yards south of the East-West fence’’, you 
will be understood. If you say: ‘‘My abscissa is 5 and my ordinate is 
-3’’, you will be using mathematical language and your less mathematica] 
friends may regard you with the uplifted eyebrow. (If you continue the 
study of mathematics, you will soon discover that mathematicians tend to 
use ‘‘trade’’ names in order to conceal the simplicity of their concepts 
and impress the uninitiated.) We now change the lettering to agree with 
common usage (Figure 2). WE becomes X'X and is called the x-axis; SV 
becomes YY and is called the y-axis. Furthermore any point P is denoted 
by P(x,y) where x, the abscissa, is the distance measured from Y'Y 
positive if in the direction of OX, negative if in the direction of OX’; 
and where y, the ordinate, is the distance measured from X'X, positive 
if in the OY direction, negative if in the OY’ direction. x and y are 
called the coordinates of point P. Then A has coordinates (5,-3); 
B, (5,1); and O, (0,0), called the origin. Thus we have made use of the 
real number system of algebra to set up a relation between algebra and 
geometry, the point being a geometrical concept. 

To Rene DesCartes (1596-1650), a famous French philosopher and math- 
ematician, is generally credited the study of plane geometry by means of 
algebra—namely plane analytic geometry. This invention also marks the 
beginning of modern mathematics. 

To illustrate the method we now formlate two geometrical problems in 
terms of coordinates.-The first proposes to find the distance between 
two points in a plane. Axes and a unit are chosen so that the coordi- 
nates of A are (3,2) and B, (7,5). The coordinates of C are (7,2). Then 
(see Figure 3) 





AC = MC ~MA= 7 - 3 = 4, 
CB = NB -N = 5 - 2 = 3, 
and 
AB = JAC’ + CB = 14? + 3% = 5, 


using the Pythagorean theorem. 
To obtain a general formla for distances, we use letters and obtain 
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(see Figure 4) 





AB = \(u - x)* + (v - y)* 
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Figure 3 Figure 4 
The second problem involves the direction of a line. This can be 


given by the size of the angle 6 (measured counter-clockwise) which 


the line makes with the positive direction of the x-axis. Hogever we 


shall use the slope m of AB which is defined as the ratio ra Hence 
AC 


- y : . 
=——— , Such ratios are in common use in measuring railroad grades, 


highway grades, slant of a roof and so on. 

The power of the two formulas just developed is wel] illuscrated by 
the following problem with which some of you may have struggled: 

Two ladders of lengths 3 yards and 4 yards (Figure 5) are placed 
against two buildings on opposite sides of an alley as shown in the 
diagram. If these ladders cross at a point 1 yard above the pavement, 
how wide is the alley? 

We choose the side of one building as OY and the pavement as OX. Then 
we assign coordinates to all the points involved, using letters to 
indicate the unknown values. Now we employ the method of algebra, 
namely, we seek to set up as many equations as we have unknowns. We 
have OB = 4 and AC = 3. Using the distance formula gives 


V(x - 0)? + (y - 0)? = 4 and V(x - 0)? + (0 — ¢)* = 3. 





slope ( 


\lso we know that the slope of OM = the slope of MB; the 
AM = the slope of M’. The use of the slope formla gives 


? 


s-O.y-1.1-0.t-1 
t= @ = S-X Q-s 


Thus we have expressed our geometric problem as one in algebra hecause 
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we have now to solve a set of four equations. Since we wish x, we use 
algebra to eliminate y, s and t and then to solve the resulting equation 
in x. If you are versed in algebra you can readily finish this problem, 

For a further illustration we return to the graph with which we 
started. Assume that the graph is a straight line and let P(x,y) he any 
point on this line (Figure 6). Then the slope of AP = the slope of AB 


Als 
its 


6(#,¥) 





c (0,t) 


mM (4,/) A(0,32) 








0 (a0) A (x,0) 











Figure 5 Figure 6 


no matter where P is as long as it is on the line. Hence 


y- 32. 86 ~ 32 , or y= 9 x + 32. 
x- 0 30 - 0 5 





of 
wht 
This means that we have found an equation which corresponds to the graph inc 
in such a way that those and only those pairs of numbers which satisfy 
the equation are the coordinates of points on the graph. We are now in 
a position to give the value of y for each value of x and vice versa. 


If x = 10.5, y = 2(10.5) + 32 = 50.9; if y = 0, x = -1%; and so on. 


and 


whe 


To express briefly the above correspondence, we say that y = 2x + 32 is the fiy 
equation of line AR, fre 
par 
OF 


uni 


In the above we can think of the straight line AB as the path traced 
by P if it moves always in the same direction. Suppose now that P, 
instead of moving along AB, moves so that it is always 5 units distant 
from C(1,2) (see Figure 7). Then the path is a circle. Using the 
distance formula we get the equation of this circle, namely: 





Thi 


pal 


V(x ~ 1)? + (y - 2)* = 5, or (x - 1)? + (y - 2)? = 25 


whence 


y=2+ J25 - (x - 1)” 
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Also (x - h)* + (y - k)® = r* would be the equation of any circle having 
its center at (h,k) and radius r. 





7 














Figure 7 Figure & 


For both the above straight line and circle, we have obtained a 
relation between x and y which is such that when x is given, y can be 
found. In such a relation x is called the independent variable, y the 
dependent variable; and y is said to be a function of x, written 
y = f(x). We can now state in general terms a fundamental principle of 
analytic geometry, namely: 

y = f(x) is the equation of a certain curve, if and only if the pairs 
of numbers which satisfy the equation are the coordinates of points 
which lie on the curve. Note that the word ‘‘curve’’ as used here 
includes a straight line. On account of the close relation between curve 
and equation, it 1s customary to abbreviate the expression ‘‘the curve 
whose equation is y = f(x)’’ to ‘‘the curve y = fix)’. Thus we say 
‘‘the line y = 9 + 32’’, ‘‘the circle x“ + y~ = 16’’ and so on. 


Suppose now that the point P(x,y) moves so that its distance from a 
fixed line MN (called the directrix) is always equal to its distance 
from a fixed point F (called the focus). Such a curve is called a 
parabola (Figure 8). To find an equation, we take MN as the y-axis and 
OF as the x-axis, OF being perpendicular to MV. Choose OF equal to one 


unit. Then for any point P(x,y) on the parabola we havd DP = FP, whence 





<= 1(x — 1)* + (y - 0)? , or y* = 2x - 1. 


This is only one of the many equations which could be obtained for this 
parabola by various choices of axes and a unit of measure. However all 
these equations would be of the second degree. 

We note in passing some applications of the parabola. The path of a 
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projectile if air resistance is neglected is a parabola. Also the ney 
television reflector on Mount Wilson, Calif., has a shape obtained by 
rotating a parabola about its axis of symmetry. 

[t now appears that al] circles and all parabolas have second degree 
equations. It can be shown that any second degree equation—namely 
Ax* + Bxy + Cy? + Dx + Fy + E = 0, where the coefficients are real 
numbers—is one of the following: circle, parabola, ellipse, hyperbol, 
or a pair of straight lines. Examples of the latter three are, respec- 
tively: 4x° + Sy? = 36, 4x* - 9y* = 36, and 4x2 - 9y? = 0, with graphs 


as shown (Figures 9, 10, 11). 














Figure 9 Figure 10 
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Figure 1] Figure 12 


It seems then that some curves can he classified according to the 
degrees of their equations. We have noted the second degree equations 
above. Also, any first degree equation is a straight line; a third 
degree equation was classified into 78 types by Sir Isaac Newton. Such 
equations as these are called algebraic curves and have an extensive 
literature. 

We have seen how algebra aids in the study of geometry. Conversely, 
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geometry illustrates some of the processes of algebra. For example, 
consider the problem of solving 2x - y=2 and x + 2y = 4 simltan- 
eously. This means to find all the pairs of values which will satisfy 


both equations. Since the equations represent straight Jines, any such 


€ new 
ed by 











Peree 

i pair will be the coordinates of a point which is on both lines and 
real therefore at their intersection. Because straight lines intersect at 
bo], only one point, we anticipate that the algebra will give only one pair 
pec- of values. We proceed to a further geometric interpretation of the 
‘aphs algebraic solution (see Figure 12). 






2x -y=2 line (a). 
x+2y= 4 line (b). 










Eliminate y and get 








x = 8 line (c). 
5 
Eliminate x and get 
= 6 line (d). 


5 
i 


Thus the solution x ==, y = 
c 






corresponds to the point [8.5 which is 


WIN 


the intersection not only of the given lines, hut also of two lines 
parallel to the y- and x-axes respectively. 

In general the algebraic solution of two equations yields the coordi- 
nates of the points of intersection of the corresponding curves. 

We shall now consider a curve which is not algebraic and which was 
defined and studied by the ancient Greeks without the aid of algebra. 
This curve, called the cycloid, is the path traced by a point on the rim 
of a rolling wheel. The figure shows the initial position and a subse- 
quent position of the wheel as it rolls along the level ground (taken 
as the x-axis). The point on the rim is initially at O and P is any 
subsequent position. 6 is the angle through which a spoke has turned and 
OP is an arc of the cycloid (Figure 13). Evidently the curve will 
consist of a series of arches each of height 2 units, the radius of the 
whee] being taken as the unit. 
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Using some trigonometry we get the equations 




























[ 
and 
An & 
call 
For each value of 6 there are values of x and y which locate a point. AOP 
Hence these two equations serve the same purpose for which we have I 
formerly used one equation. They are called parametric equations of the by F 
curve and 9 is called a ‘‘parameter’’. To obtain one equation we elimi- 
nate @ and get lin 
sitalilancininecigy num 
x = arc cos (1 - y) —V2y - y* i 
However the parametric equations are the more useful. With the aid of cur 
the calculus they enable us to study vehicular motion—motion of an par 
automobile for example. sun 
If the arch of the cycloid is inverted, then three particles sliding 
without friction from O, A and B at the same time will arrive at B, the 
lowest point, simultaneously (see Figure 14). Also a particle sliding 
without friction will go from O to B in less time than along any other 
curvd joining O and B, 
The coordinates used thus far are often called rectangular because, 
if lines are drawn through the points parallel to the axes, a network of 
rectangles is formed. However a point in a plane can be located by means 
of to numbers in mary other ways. We mention one of these ways. 
A line Q’@ is drawn through a fixed point O (called the pole) making 
an angle @ with a fixed line OM (called the polar axis). A distance r 1s 
measured from O along OQ to point P. in case r is positive (Figure 15). 
If r is negative, it is measured along 0Q' to P . Thus a point P is 
na 
p(r,@) " 
na 
a 
ha 
of 
- 
Mm 
Figure 16 
Bi 
located and P is said to have the polar coordinates r and 0, written 7 
P(r,@). For example (3,30) is at A; (-2,30°) at B; (3,0°) at C and 
so on. 
Just as y = f(x) represents a curve, so also does r = f(@). As an r 


example consider the curve defined as follows: 
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Draw a circle of radius | and center at (1,0). Now draw any line OQ 
and extend it 1 unit to P(r,@). Then OP = 00+ @, or r = 2 cos 6 + 1. 
An arc APS of the curve is shown (Figure 16).This curve is sometimes 
called the trisectrix because, no matter where P is on the curve, angle 
AOP is one third of angle MOP. 


For another example we shall find the equation of the curve traced 


by P(r,@) if = = e, where DP is the perpendicular distance from a fixed 


line AB a units from O and perpendicular to OM, and e is a positive 
number less than 1. Then OP = e « DP, whence r = e(a +r cos @), and 


For various values of a and e we get a group of 


curves, but it turns out that they are all ellipses (Figure 17). In 
particular the path of the earth around the sun is an ellipse with the 


sun at O. 














Figure 17 Figure 18 


As we have seen, an ellipse has an equation in rectangular coordi- 
nates, but for the purposes of astronomy it turns out that the polar 
coordinate equation is superior. This raises the question ‘ ‘What coordi- 
nate system is best?’’ Since there are many systems,probably there is 
a hest one for each geometrical problem. However the two systems we 
have mentioned are the simplest and are applicable to a large number 
of problems. 

It is interesting to compare the polar coordinate equation of a curve 
with its equation in rectangular coordinates. To do this we let the 
pole and origin coincide, alsothe polar axis and the x-axis (Figure 18). 


is - : 2 2 2 . i 
y inspection we have r = a anc x + y =a ~= as equations of the circle 
with center at O and radius a. Evidently r = a is the simpler. 


On the other hand the straight line y= : 39 becomes 





é ,?” . 
—— 160 , a much more ‘ ‘unhandy form. To make this 
2 sin G - 9 cos @ 
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‘*transformation’’ we note from the figure that x = r cos @ ang 
y = r sin &, substitute in the rectangular equation, and do the 
necessary algebra. 

In conclusion it should be said that a person familiar with plane 
analytic geometry may say that this discourse is noteworthy (if at all) 


for what has been left ovt. The reader may decide this question by 


studying a text on the subject or by taking a course in it. In either 


Je to fill in the algebraic and trigonometric details 


which the author has so considerately omitted. 


event he will be al 
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THE STRAIGHT LINE TREATED BY TRANSLATION 
AND ROTATION 


Kenneth May 


Transformations are usually introduced rather late in the first 
course in analytical geometry and appear as hardly more than a device 
to facilitate the graphing of conics. What should be a fundamental 
notion is merely a footnote. The following remarks indicate one way in 
which an early introduction of the equations of translation and rotation 
can be turned to immediate account in developing the basic theorems on 
the straight line. At the same time it appears that the usual proofs 
based on the geometrically established constancy of the slope are 
redundant. The nature of the graph of Ax + By = C follows from the 
linearity of the axes without reference to the notion of slope, which 
can then be introduced and proved constant analytically. 

Consider the equation Ax + By = C, where if A= 0, B# 0; if B= 0, 
A# 0; or A¥# 0, BA 0. Suppose A = 0. Then if C = 0, the graph is 
evidently precisely the x-axis since there and only there is y = 0. If 
C= 0, we translate the origin to the point (0, 5), and the graph is now 
seen to coincide with the x-axis. The case B= 0 is treated similarly. 
If neither A nor B is zero, we translate the origin to the point (©, 9), 


and the equation becomes Ax + By = 0. Rotation through an angle @ yields 


(AcosO + Bsin@)x + (Bcos@ - Asin@)y = 0. 


Evidently the first term vanishes if tan@ = - . Moreover, for this @ 


the coefficient of y cannot vanish. Hence the equation becomes y = 0 
and its graph is coincident with the x-axis. 

Thus Ax + By = C is shown to be the equation of a straight line since 
it is in all cases congruent to an x-axis. Moreover, the argument 
indicates the relation between the coefficients of the equation and the 
slope and inclination of the line. The special forms of the equation and 
(y2 - y;) 


follow immediately. 
(x, - x,) 


the constancy of the ratio 
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A course in intermediate algebra designed for 
college students who did not study a second 
course in algebra in high school. Appropriate 
refresher work in arithmetic is included and 
throughout the text there is continual emphasis 
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Comment On Two Papers By P, A, Piza. By H. W. Becker 








I wish to congratulate Senor Piza upon his elegant theorems in 
escalator number theory, and to flash some sidelights. 

Where {2 and 22 are any two operations, s'lz = s{0'z is an operational 
equation. The escalator numbers A, are then solutions of an op. eq. 
between + and X (1). These with exponenting are the first three of an 
infinite series of operations, which have been variously defined (2-6). 
These operations each have an inverse (as / is the inverse of X) and a 
dual (as || is the dual of +, thus s||z = sz/(s + z) is the dual of 
s + z). One innovator (7) proposed that we write the third operation 












. ° ° ° 1 fa 
s* = s$z, because ‘‘money is power’’; and its inverse s /z = s¢z, 


because ‘‘money is the root of al] evil’’. Howsoever that may be, + = $ 

would he a concise instruction to solve an op. eq. s + z = s§$z, etc. 
The fundamental equation of escalator number theory is "a © Ags 
x et ee = + = ,? ° = 

A A. , whence 1=AA /(A, Rat AIIA,4, . Or in elec 


n 1 non+i 
trical parlance, ‘‘The joint impedance of two successive cobasic 














escalator numbers in parallel is an impedance match of unity 

The impulse behind this theorem goes back to the very dawn of Number 
Theory. That 1 +2 + 3 = 1 x2 x3 so impressed the ancients, they called 
6 the first of the ‘perfect numbers’. They went so far as to inscribe 
the related ‘amicable numbers’ on pellets, and eat them for an aphro- 
disiac (8). 

Being with Piza a connoisseur of recurrence gems, I wish to exhibit 
further references for his excellent development of the Kummer Numbers— 
called cummlative numbers, and similarly employed (9). They perhaps 
Originated in a summation of series by Euler (10) in a letter to 

| Stirling, dated ‘'27 Julii 1738’’. 
They have a combinatory interpretation (11), :K. is the number of 
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permutations of t letters with c - 1 inversions. E.g., abc has uy 


inversions, cba has two, and bac has one inversion of standard orde; 


between consecutive elements. They also enumerate ascending runs anj 
triangular permutations (12). They have yet a different type of inter. 
pretation: ,K_ is the number of words of t letters, c different, such 
that the pth position may be occupied by any of the first p letters of 
the alphabet. Thus aaa is a member of ,K,, aac of ,K,, and abc of 4k. 


The expositor’s enthusiasm for the summational advantages of the] 
formula is shared by practical statisticians (13-16), among whom it has 
largely superseded Bernouilli’s method (17)—which discourse wil] hoy. 
ever remain long a classic in the literature of jubilation, and symbolic 
polynomials. 

Devotees of both music and mathematics may be interested in the 
combinatory interpretation of ,D. = A‘o*, the number of musical cadences 
of t notes, c different. Thus that most celebrated example bbba, the 
opening bar of Beethoven’s Fifth Symphony, is a member of ge = 14. For 
t < 6, these may be written down by permuting all the rhyme schemes of 


page 23, XXII, this Magazine, since -. =c! 7. 
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Tables of Bessel Functions of Fractional Order, Vol. I., Columb 
University Press, N. Y., 1948. 413 pp. Price $7.50. The Computation 
Laboratory, National Applied Mathematics Laboratories, Nationél 
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Rureau of Standards. 


This is another of the series of tables produced by what was orig- 
inally known as the NYMTP (New York Mathematical Tables Project) which 
hecan under the auspices of the Works Progress Administration for the 
City of New York. The first table appeared in 1939. Work continued with 
the aid of funds from the Office of Scientific Research and Development 
and was operated by the National Bureau of Standards. The project has 
now become one of the National Applied Mathematics Laboratories. However, 
it is to miss the whole truth if one looks only at the agency from which 
the table has emerged and forgets the human side of the picture, the 
men behind the table: They are Arnold N. Lowan, now chief of the labora- 
tory, and M. Abramowitz, G. Blanch, A. Hillman, W. Horenstein, M. Karlin, 
J. Laderman, I. Rhodes, H. E. Salzer, I. Stegun of the technical staff. 
It is out of their drive, hopes and feelings that tables and programs 
for tables arise and have their being. 

In the growing field of mathematical technology, mathematical tables 
and other aids to computation assume central importance. Only those 
functions that are usable without tabulation or that are adequately 
tabulated can be used in the solution of practical problems, requiring 
numerical answers. Numerical results are the basis for comparison of 
theory and measurement and for the prediction of the performance of a 
proposed device. When a function has been tabulated, e.g., sin x, cos x, 
log x, etc., it becomes available for the expression of the solution 

a mathematical-technical problem in ‘‘closed’’ form. Almost al] 


problems in mathematical technology which involve the deduction of the 


consequences of laws of physics, also involve approximations and 
simplifying assumptions. In many cases, the factor which governs these 
simplifying assumptions is the existence of known or tabulated functions 
to which the assumption leads. Even the use of Fourier series and power 
series can be viewed as an approximation which yields the answer as a 
linear combination of trigonometric functions or as a polynomial. 

The foreword by Schelkunoff, the introduction by Abramowitz and the 
bibliography are important parts of the table. The foreword indicates the 
lifferential equations which are solvable in terms of the Besse] 
functions of fractional order. It is also pointed out that among the 
problems which lead to these differential equations are wave propogation 
in stratified media with a constant gradient in the index of refraction 

a non-uniform transmission line. The introduction discusses the 
mathematical properties of the functions tabulated, differentia] 
equations reducible to Bessel’s equation, zeros, interpolation, Besse] 
functions considered as functions of their order (by Salzer) , method 
of computation, and the preparation of the manuscript. The bibliography 
flves thirty references: a) tables and texts, 1-ll, 6) articles relating 
0 the theory of Pessel functions, 12-22, c) texts and articles relating 


to physical applications, 23-30. 
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The technical details of the table are: 
1. Tables of J, (x) 
Vv x 


- 3/4, - 2/3 [0(.001)0.9(.01)25; 10D) 

~ ifs, - 1/4 [0(.001)0.8(.01)25; 10D) 
[ J 
[ 





1/4, 1/3 0(.001)0.6(.01)25; 10D 
aie, 3/4 0(.001)0.5(.01)25; 10D) 





2. Tables of A,(x) and B,(x): 


J,(x) = A,(x)cos(x - m/2 - 77/4) + By (x)sin(x - m/2 - 7/4) 


V Vv 
= 1/4, 1/3, 2/3, 3/4 
x = [25(.1)50(1)500(10)5,000(100)10,000(200)30,000; 10D] 


3. Zeros of J (x) 


4. Tables of L,(u)—For interpolating in the V-direction. 


5. Other auxiliary tables. 


There is no question of the reviewer judging the merits of the 
tables; they are excellent tables and fill an important need. Those who 
are interested in the general subject of tables will find two sources 
helpful. The first source is the journal ‘‘Mathematical Tables and other 
Aids to Computation.’’ The second source is the book ‘‘An Index of 
Mathematical Tables,’’ by A. Fletcher, J.C.P. Miller, and L. Rosenhead, 
publisher McGraw-Hill, N.Y., 1946. 


Concord, Massachusetts Nathan Grier Parke II 





Fundamentals of Statistics. By J. B. Scarborough and R. W. Wagner. 
Boston, Ginn and Co., 1948 7 plus 145 pages. $2.40 


Books in statistics written by statisticians usually turn out well 
while those written by non-statisticians usually turn out badly. 
Unfortunately practically all elementary books in statistics belong in 
the latter category; the present volume was written by two mathematicians 
at the U. S. Naval Academy. Apparently they made the mistake of letting 
themselves be guided by other elementary books in the field. 

The book is intended as a textbook for a brief course in statistics 
for students who have had calculus. It is clearly written and contains 
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a wealth of exercises. Considerable material is packed into its 108 
pages: frequency distributions, the first four moments, correlation, 
regression, binomial and normal distributions, Gram-Charlier series, 
and a few large sample tests of significance. Also there is an appendix 
on combinatorial probability, and tables of the normal distribution and 
and its first four derivatives. 

The trouble with the book is in its selection of materia]. The main 
problem of statistics is uncertain inference: sampling theory, esti- 
mation theory, and tests of hypotheses. The book relegates all these 
fundamental problems to a brief final chapter. It chooses instead to 
emphasize the frequency tally, the class interval, the miscellaneous 
means, skewness and kurtosis, computational details, and the Gram- 
Charlier series. These topics are included at the expense of Student’s 
distribution, the Chi-square distribution with al] its many uses, and 
the most powerful technique the statistician has—the analysis of 
variance. The student who surveys statistics through these pages will 
not see it as the exciting subject it is, but merely as an insane 
accounting system apparently bent on reducing any number of figures to 
four figures. 

Although it now appears utterly useless to do so we are duty-bound 
to observe some of the ancient errors reproduced in this book. First of 
all there is the nonsense about skewness and kurtosis. Statisticians 
have been poking fun at skewness and kurtosis in their classrooms for 
years. They have buried skewness and kurtosis in their journals time 
and time again. Yet here they are again, not only hale and hearty, but 
of sufficient stature to be included in a severely limited selection of 
topics. 

Another thing statisticians like to exercise their sarcasm on is the 
band about a regression line bounded by parallel lines. This text does 
not deprive them of that opportunity. 

A favorite classroom sermon has to do with the sin of equating 
parameters to their estimates. The authors do just that on page 105. To 
me, page 105 is the saddest part of the book for here the authors commit 
two mathematical crimes. The second one is in leading the student to 
believe that an unbiased estimate of a function of a parameter may be 
obtained by using the same function of an unbiased estimate of the 
parameter. 

Of course every book has its little errors, and by emphasizing them 
a reviewer can give a bad impression of a good book. But the errors 
mentioned here are not of that category. Statisticians will immediately 
deduce that authors who can make such errors do not read the statistical 
journals, have not read Fisher’s books, have not read the Neyman and 
E. S. Pearson papers, have probably not even read Kar] Pearson. In short 
they have attempted to write a textbook in a field about which they 


are 1]]-informed. 


A. M. Mood 
The authors will reply in the May-June issue of this Magazine—fd. 












PROBLEMS AND QUESTIONS 
Edited by 





C. G. Jaeger, H. J. Hamilton and Elmer Tolsted 


This department will submit to its readers, for solution, problems which seen 
to be new, and subject-matter questions of all sorts for readers to answer or 
discuss, questions that may arise in study, research or in extra-academic 
applications. 

Contributions will be published with or without the proposer’s signature, 
according to the author’s instructions. 

Although no solutions or answers will normally be published with the 
offerings, they should be sent to the editors when known. 

Send all proposals for this department to the Department of Mathematics, 
Pomona College, Claremont, California. Contributions must be typed and figures 
drawn in india ink. 


SOLUTIONS 


No. 21. Proposed by Julius Sumner Miller, Dillard University 


A plane is inclined to the horizontal at an angle B. At the foot of 
this plane a particle is projected with velocity V at an angle A with 
the plane. Find the condition for maximum range. 


Solution by Howard Eves, Oregon State College. 


This is an old problem; see, e.g., J. H. Jeans, An Elementary 
Treatise on Theoretical Mechanics, pp. 209-210. It is there shown that 


o §& 
A= 45° -3 


That is, to get maximum range up the inclined plane we must project in 
the direction which bisects the angle between the inclined plane and 
the vertical. This maximum range has the value 


y2 
(1 + sin B). 





Solved also by Leo Moser, Winnepeg, Canada. 


No. 27. Proposed by C. . Trigg, Los Angeles City College. 


In a square of side “a” a right triangle is inscribed with one acute 
angle at a corner of the square and the other two vertices on the sides 
of the square non-adjacent to that corner. If one of the legs of the 
triangle is equal to the other plus the side of the square, find the 
longer leg of the triangle. 


Solution by C. D. Smith, University, Alabama 


I first solve the more general problem for a rectangle with sides 4 
and b. (a > 6). The square is a special] case. The right triangle has 


vertex of the right angle along side b with legs c and d. (c > d). 
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Let x be the third side of the triangle formed by c and a. The angle 
formed by c On a is equal to that formed by d on b, hence the relation, 


cla = d/(b - x). Also x? = c* ~ a*, Eliminate x and apply the given 


condition d = c — a. Using the relation r = b/a we have by obvious 






reductions the equation 








(1) os ee tt s op vy + Bee 0? = 7) e « &. 






A positive value of c greater than a is a solution. For r = 1 we have 


the special solution for a square. 





—— 


a -_ 
The positive root of (2) is c “20 + V5). 


























Jalved also by the proposer, and Adrian Struyk, Patterson N. J.; 4. B. 
Clarke, San Jose; and Francis L. Miksa, Aurora, Illinois. 


PROPOSALS 


34. Proposed by Victor Thebault, Tennie, Sarthe, France. 









Find a number of four digits such that its square ends with the same 
four digits in the same order. Show that its cube and its fourth power 
have the same property. 








35. Proposed by Victor Thebault, Tennie, Sarthe, France. 





A plane P divides the volume of a sphere into two parts V and V’ and 
determines two spherical segments of areas S and S’. If it is known 
that S/S’ = k, calculate the ratio V/V’ in terms of k. 











36. Proposed by Julius Sumner Miller, Michigan College of Mining and 
Technology. 





A gun can put a projectile to a height R/n, R being the radius of 
the Earth. Assuming the variation in gravitational force with altitude, 
find the area commanded by the gun. 












Problems 18 and 19 in Vol. 22, No. 1 were proposed by Julius Sumner 
Miller rather than by Julius Sumner as there indicated. Problem 18 as 


solved in Vol. XXII No. 2 should have been credited to K. L. Cappel. 
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Edited by 
Marion E. Stark 


Let us know (briefly) of unusual and successful programs put on by your 
Mathematics Club, of new uses of mathematics, of famous problems solved, and go 
on. Brief letters concerning the MATHEMATICS MAGAZINE or concerning other 
“matters mathematical” will be welcome. Address: MARION E. STARK, Wellesley 
College, Wellesley, 81, Mass. . 





NBS Interim Computer 

“Design and development work supported by the Air Force is wel] under 
way at the National Bureau of Standards for the construction of a small]- 
scale electronic computing machine to be used until the several! large- 
scale machines now being built become available. The new high-speed 
machine, to be known as the NBS Interim Computer, will perform as sub- 
stantial portion of the computation work of the Bureau’s laboratories, 
solving many problems until recently considered impossible of solution. 
It will also aid in computing machine development at the Bureau and will 


provide important training and operational experience for personnel of 


those agencies that plan to operate the more complex electronic con- 
puter as soon as their construction is ¢ omplet e 
seem : ; 
The National Bureau of Standards is now engaged in an extensive 


computer program, undertaken in cooperation with the Office of Naval 
Research, the Bureau of the Cen us, the Department of the Army, and the 


‘ 


Department of the Air Force. This program involves 


the research, design 
and development work necessary to produce electronic machines that wil 
perform, upon instruction, predetermined sequences of calculation 
running into the thousands of operations without the intervention of 
human operators The result will be the solution in a few hours of 


complex problems in atomic physics, ballistics, and aero-dynamics which 
cannot now be solved except by simplifying assumptions and thousands ol 


man-days of work. The rapidity with which numerical data can be handled, 

. . . - ” 
classified, and analyzed will also be correspondingly increased. 
National Bureau of Standards Technical Report 1321] 
Washington 25, D. C. 





Sums of Two Squares 
; 2 2 
Beginners in algebra are taught the factors of a - 6 , but are 
usually impressed with the fact that the sum of two squares cannot be 
factored. Certainly a* + b* has no real factors of the type considered 


. 9 
in elementary algebra. But (a*)* + (2b7)? has the factors 


2 ¢ 2 ‘ ‘ 
(a2 + 2ab + 2b7)(a2 - 2ab + 267). 
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There are many cases in which a numerical integer, found as the sum of 
two squared integers, has integral factors. The search for such integers 
raises some interesting questions. 
The we!l known formula which gives integer sides for right triangles, 
(a2 - b2)? + (2ab)2 = (a2 + b?)(a? + 67) 
gives two equal factors for a sum of two squares. 
If we ask what integers are the sum of two integers squared in two 


} " - 2 
different ways, we start with a® + b* = x? + , whence a‘-x* = y*-b ; 


yeas cd? + OH. 0 


This equation can be satisfied by taking the product as 4PQand making 


a + x = 2PQ y + b = 2P 

a-x=2 y - 6 = 2Q 
whence 

a = PQ + | b=P-Q 

x = PQ-1 y=P+Q 


So we have two cases of the sum of two squares and the factors 
(PQ + 1)? + (P FQ)? = P2Q? + Pe + G+ 1 = (P? + :1)(0* + 1) 

For example, if P = 2, Q = 3, this formula gives 

Pete «+ Ss Sx 10° 


‘ 


50 is the smallest number that can be expressed as the sum of two 


squares in two different ways. 
To make sure that a and b differ from x and y, note that a # x, and 


ifa# y, PQ+ 1A P+ Q, which gives (P- 1)(Q-1) 4 0, so it is only 
necessary that neither P nor Q equals 1. 


The formla P?Q* + P? + Q + 1 may be taken as the sum of two squares 
in a third way, if we impose an additional condition. This formla ends 


in a square, 1, and P*Q* + P* + Q* is a square if (P*)* = 4(P2Q* )Q?, 
that is if P = 2Q*. Using this value we get 

4h + 49% + Q? + 1 = (2Q% + Q)? + 17 = (4Q% + 1)(Q* + 1) 
This gives us another case of factors for the sum of two squares. 


This product, (4Q* + 1)(Q* + 1), provided Q# 1, gives a number 
which can be separated into the sum of two squares in three different 


ways: the original two ways are 


(203 + 1)2 + (292 FQ)? 


and the third way is 
(2Q° + Q)? + 1? 


For example, if we take Q = 2, we get 
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(16 + 1)? + (8 = 2)? = (16 + 2)2 + 1% = 65 x 5 = 325 


325 is the smallest integer that can be separated into the sum of two 
squares in three different ways. 


Tufts College William R. Ransom 





More news from France, sent us by Colonel Byrne. Professor Gaston 
Julia has been elected Vice President for 1949 of the Academie des 
Sciences. Three volumes of the collected works of Poincare have been 
published, while seven more await publication. Inflation in France has 
reached such a state that twenty million francs will be needed for that 


publication. 





New York, January 27: A feature of the annual meeting of the Duodecimal 
Society of America, held tonight at the Gramercy Park Hote] in New York 
City, was the explanation of a new, mathematically—hased notation for 
music by Velizar Godjevatz. Mr. Godjevatz, who has studied music at 
Oxford, Grenoble, and in Berlin, pointed out the illogic of applying 
to our present 12-tone tempered system the staff of five lines and the 


concepts of the octave and a 7-tone system. 

“‘It should be borne in mind,” he said, “that our present musical 
notation was conceived for another musical system in use centuries ago. 
When Werckmeister definitely established our 12-tone equidistant system, 
and Bach drew the necessary conclusions in publishing his Well Tempered 
Clavichord in 1722, the notation was not changed to accomodate the 
practical results obtained by embracing this new system.” 

He proposed the “Notation Godjevatz” in which each of the twelve 
tones within a Region (his term for the sound range of an octave) has 
seperate representation on a staff of seven lines, with the middle line 
eliminated for easier scanning except when a note falls upon it. This 
staff requires neither sharps nor flats, and is proposed as greatly 
simplifying the learning and the performance of music, both for the 
amateur and the professional. 

Mr. Godjevatz explained the mathematical relationships among musical 
tones, calling a notation that would make these relationships plain the 
most important problem in music today. “It would be advisable to employ 
for our dodecaphonic musical system the duodecimal number system,” he 
declared. ‘‘Every tone could be represented by one digit only, and this 
would be its relative pitch. The adoption of the duodecimal musical 
notation would create a sort of msical arithmetic. This would encourage 
more people to learn to play an instrument, people who otherwise would 
not expend the time and effort required to master a good reading 


knowledge of music.” 








